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Abstract 

We study the effects of a uniform magnetic field on the one-dimensional spin- 
orbital model in terms of effective field theories. Two regions are examined: 
one around the SU(4) point (J=K/4) and the other with K <C J [see Eq. 
(1)]. We found that when J < K/4, the spin and orbital correlation functions 
exhibit power-law decay with non-universal exponents. In the region with 
J > K/4, the excitation spectrum has a gap. When the magnetic field is 
beyond some critical value, a quantum phase transition occurs. However, the 
correlation functions around the SU(4) point and the region with K <C J 
exhibit distinct behavior. This results from different structures of excitation 
spectra in both regime. 
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I. INTRODUCTION 



The interest in studying the role of orbital degrees of freedom stems from the understand- 
ing of the magnetic structures of transition metal compounds [f[l§. In these systems, the 
low-lying electron states have orbital degeneracy as well as the usual spin degeneracy. This 
may result in interseting magnetic properties of the Mott insulating phase. For example, 
the magnetic ordering is influenced by the orbital structure which may change under the 
pressure or the magnetization is a nonlinear function of the magnetic field even in the case 
of an isotropic exchange interaction of the type Si ■ Sj |j| . These distinguish the spin-orbital 
models from the ordinary Heisenberg model with the spin degrees of freedom only. Thus, to 
understand the magnetic properties of these compounds, an investigation of the interplay 
between spin and orbital fluctuations is necessary. 

As a prototypical model in which the quantum fluctuations of the orbital degrees of 
freedom are important, we would like to consider the following Hamiltonian: 

H = Y\K{Si ■ S l+1 )(f t ■ f l+1 ) + Jj t ■ S l+1 + J 2 f, • f i+1 - hS?], (1) 

i 

where Tj is the pseudospin to represent the orbital degrees of freedom. We assume that K 
and h are positive. The Hamiltonian ([!]) is related to the recently discovered spin-gapped 
materials, Na 2 Ti 2 Sb 2 and Na 2 V 2 05 [|J. These materials have a quasi-lD structure 
and are modeled by the quater-filled two-band Hubbard model which is equivalent to Eq. 
(|l|) (h = 0) in the strong Coulomb repulsion limit. Instead of a magnetic field, we can also 
impose a uniaxial pressure P which gives rise to a term of the type APT? in the Hamiltonian. 
In other words, a uniaxial pressure affects the orbital sector in much the same way that a 
magnetic field affects the spin sector. Therefore, our results also apply to this case where 
the roles of spin and orbital operators interchange. 

The Hamiltonian (JI|) with h = is invariant under independent SU(2) rotations in the 
spin and orbital spaces. When J\ = J 2 = J, it is also invariant under the exchange between 
S and T. We restrict our consideration to this more symmetric case and comment on the 
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asymmetric case in the conclusion. The phase diagram has been studied in Ref. [0,0] and the 



properties of the ground state depend on the ratio J/K. When J > K/A, the ground state 
is doubly degenerate with alternating spin and orbital singlets. The excitation spectrum has 
a finite gap. In the region with —K/A < J < K/A, it is a critical theory. The point with 
J = K/A is special because the symmetry is enlarged to SU(4) [^,0 and it is Bethe-ansatz 
solvable 0,0]. There are three gapless bosonic modes with the same velocity. This implies 
that the central charge c = 3 and as shown by Affleck ||15| , the fixed point Hamiltonian is 
described by the SU(4)i Wess-Zumino-Novikov-Witten (WZNW) model. The predictions 
of the conformal field theory have been confirmed by numerical work |]||]. Away from the 
SU(4) point, the low energy excitations of the Hamiltonian ([!]) are described by different 
effective theories at different range of J/K. As shown in Ref. 0, an appropriate starting 



point at weak coupling regime (J 3> K) is two decoupled Heisenberg chains [|16j while near 
the SU(4) point, we have to begin with the 0(6) Gross-Neveu (GN) model. 

In the present paper, we study Eq. ([]]) on the symmetric line with h ^ 0. The case with 
J = K/A has been done in Ref. |]TT[ with Bethe ansatz and numerical methods. Although 
the magnetic field is not directly coupled to the orbital degrees of freedom, it still affects the 
orbital structure. In one dimension, this is reflected on the change of the orbital correlation 
functions. For J < K/A, the excitations are gapless at any finite h and the leading behavior 
of the correlation functions is the same as that at J = K/A. The orbital correlators is hardly 
influenced by the magnetic field except the 2k p component. Not only the characteristic 
momentum becomes incommensurate but the corresponding exponent is nonuniversal and 
dependent of the magnetization now. Especially, it shows a discontinuous jump compared 
with that at vanishing magnetic field. This is because a marginally irrelevant operator at 
h = becomes a marginal one under the magnetic field. 

For a gapped spin liquid, there is a quantum phase transition induced by the magnetic 
field. When h > h c (the critical field), the magnetization M ^ whereas M = for 
h < h c . The corresponding quantum critical point (QCP) is determined by the gap of 
the S z = ±1 components. When the magnetic field exceeds this gap, the corresponding 



excitations become gapless due to the condensation of these S z = ±1 bosons. For J > K/A 
but still near the SU(4) point, the excitations with S z = are still gapped and we conjecture 
that they are described by the 0(4) GN model. This results in a change on the spectrum 
when M ^ 0. For h > h c , the low energy excitations with S z = are kinks while for 
h < h c , they become massive fermions which can be considered as the bound states of kinks. 
In this case, the orbital correlators are exponentially decaying functions with or without 
algebraically decaying prefactors while the spin correlators show algebraic decay. However, 
when it comes to the region with K <C J, the behavior changes though the ground states 
in both regions have similar properties. Now the quantum phase transition occurs at the 
value of twice the gap instead of the gap. In addition, both types of correlation functions 



exhibit power-law decay with universal exponents. Indeed, as suggested in Ref. PJI4, the 
excitation spectrum at K/A < J < K/2 is different from that at J > K/2. In the former 
case, the structure factor has an incoherent background with the top at q = n and a coherent 
magnon peak at q = -k/2. The relative amplitude of the peak at q = ir/2 with respect to the 
incoherent background at q = ir decreases with increasing J/K and vanishes at J = K/2. As 
for the latter, there is only an incoherent background at q = tt |T7|]. Our analysis indicates 
that both regions have distinct responses to the magnetic field (or uniaxial pressure). 



II. THE EFFECTIVE HAMILTON! AN AROUND THE SU(4) POINT 

We start by a brief review of the derivation of the low energy effective Hamiltonian to 
fix our notation. Following Ref. 0, the low energy effective Hamiltonian of Eq. ([!]) (h = 0) 
at J = K/A can be derived by considering the following repulsive SU(4) Hubbard model 
(U > 0) at quarter-filling: 

Hu = -tJ2( C i+laa C iaa + H.C.) 
iau 

+77 J2 n iaan iba '(l -5 ab 5 aa >). (2) 

iabaa 
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Here cf acr creates an electron with the orbital index a = 1,2 and spin a =f, I and n iaa = 
cf aa c iaa . The quarter-filling electron band implies that the Fermi momentum k F = where 
a is the lattice spacing. 

The low energy physics can be described by the right-moving (R aa ) and left-moving 
(L aa ) fermions, which is related to the original lattice electron operator as: Ci ar7 /y/ao = 
Raa(x) exp {ik F x) + L aa (x) exp (—ik F x) where x = ia . We introduce bosonic fields 4>R(L) aC r, 
which satisfy the commutation relation [4>Laa,4>Rba'} = ~\^ab$ cja 'i to bosonize the right and 
left movers as: R(L) aa = (2na )~ 1 ^ 2 r] a(T exp {±.i\f^n(j) R ^ L ) aa ) where T] aa are Klein factors. It 
is more convenient to employ the following basis: 

$ c = (<M + 0ii + <M + <M)/2, 

$>s = (<M - 0i| + <M - 02.0/2, 

$/ = (01 T + 01| - 02T - 02|)/2, 

® s f = (0it - 0i| - 02T + 02|)/2. (3) 

In this new basis, $ c represents the total charge degree of freedom while other bosonic fields 
$ a (a = s, f, sf) correspond to the spin-orbital degrees of freedom. 

Substitution of the above representations into Eq. (0) gives a Hamiltonian which consists 
of decoupled charge and spin-orbital degrees of freedom. The charge sector is described by 



a Gaussian model of $ c perturbed by an Umklapp process: cos (JIGtt/ K c $ c ) where K c is 



an increasing function of U at weak coupling. As shown in Ref. ||18| , there exists a critical 
value U c (~ 2.8) such that K c = 2 where a Mott transition occurs and the system becomes 
an insulator when U > U c . In the following, we assume that we stay in the insulating phase 
and focus on the spin-orbital sector. By introducing six Majorana fermions: 



(f 1 + *£ 2 W) = -jL= exp {±iV^<S> sR{L) ), 
(£ 3 + ^ 4 W) = ^exp(±«v / 4^$ / i ?( L)), 
(e + = -p= exp (±iV^$ sfRiL) ), 



(4) 
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where r\ a (a = 1,2,3) are Klein factors JOJ, The corresponding Hamiltonian becomes 

n = --v f Y.&r R - ami) + g 3 (E w 2 , (5) 

Z a=l a=l 

where G3 = —Uao/2. This is an 0(6) GN model with marginally irrelevant 4-fermion 
interaction due to G3 < 0. Therefore, the low energy behavior of the SU(4) spin-orbital 
model is described by the SO(6)i (SU(4)i) WZNW model. 

Away from the SU(4) point, i.e. J 7^ K/4, there is an additional term in the Hamiltonian: 

AH = -(d x $ s ) 2 + G(t CrCl) 2 ~ — cos (V4^$ s ), (6) 

7T Q=3 7Ca 

where G = c(J — K/4) and c is a positive constant. With the above, the low energy effective 
Hamiltonian of Eq. ([I]) with J\ = J 2 becomes 



H = H s + H so + H, 



inti 



where 



H a = ^[(^0 S ) 2 + (d x <S> s ) 2 ] + 9 -±(d x $ s ) 2 - JLd x Q t 

I 71 \ 71 



a=3 a=S a=3 

H in t = -— cos {^ s ){g 2 J2 ZUl + 23^1), (7) 
ira a=3 

with gi = G3 + G and g% = 2G3. Here, the velocity of $ s is represented by a different 
notation. Because under the magnetic field, the spin SU(2) rotational symmetry and the 
exchange symmetry are broken. We do not expect that all fields have the same velocity. 
Note that only $ s is directly coupled to h. Physically, this reflects the fact that $ s represents 
the S z = ±1 components of magnons while other fields have S z = 0. Eq. (^) is our working 
Hamiltonian. 

Let us consider H so first. The one-loop RG equations of coupling constants are 

dg± g± 



dl 71 



(8) 
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where g± = g% ± ^2/2 and I is the logarithm of the length scale. The parameter space is 
divided into three regions where the boundaries are the same as those in the case without 
the magnetic field |J: (i) Region I (G < 0). Both g± flow to zero couplings, (ii) Region II 
(0 < G < — 2G3). g + flows to zero coupling while g_ flows to strong coupling, (iii) Region 
III (G > —2G 3 ). Both flow to strong couplings. Based on the RG flow, we discuss the 
properties of the Hamiltonian fl7|) in the following sections. 

III. CORRELATION FUNCTIONS AROUND THE SU(4) POINT 

We now discuss the spectrum of the Hamiltonian (|7|) and the behavior of correlation 
functions around the SU(4) point. 



A. J < K/A 

Eq. (H) tells us that H so describes four free massless Majorana fermions at low energy 
when J < K/4. Integrating out the high energy modes of these fermions generates two 
kinds of interactions through H int : (d x $ s ) 2 and cos (\/167rZ$ s ) where I is an integer. (It also 
renormalizes v °.) The former can be absorbed into the g\ term of H s . In addition, it consists 
of higher powers of gi and g% and we do not expect this will change the sign of g±. The latter 
(the cosine term) is associated with the lattice translation symmetry of the original lattice 
model: <E> S — > $ s + and — > — C^Cx (a, = 3, • • • , 6). They are irrelevant operators. 

Therefore, to calculate the long distance behavior of correlation functions, we can treat H s 
and H so as independent sectors and throw away H int . 

Keeping the above approximation in mind, H s can be diagonalized as the following: 

1 



H s = — 
2 



a 



id x Q s y + a (d x $ s y 



(9) 



where a = «/l + ^% and v s = av®. We also translate $ s as: $ s — > $ s + — x. With the 
definition of the magnetization M = 2Sf ot /N (0 < M < 1) where is the number of lattice 
points, we obtain M = ^«a. (This relation is valid when M <^ 1.) In general, we expect a 



and v s depend on the magnetization or magnetic field. Moreover, a < 1 because g\ < in 
this case. 

Now we can calculate the spin and orbital correlation functions. The spin and orbital 
pseudospin operators are, respectively, defined as: Si — % J2 a c ta® c ia an d = | J2a c ta^ c i^ 
where a and f are Pauli matrices in the spin and orbital spaces, respectively. Near the 
SU(4) point, they can be expressed by the WZNW fields as: Si/a = J s {x) + (e 2ikFX N s (x) + 
H.c.) + (-l)«on s (i) and Tj/a = Jt{x) + (e 2lkFX N t (x) + H.c.) + (—l)°Qn t (x) where x = ia . 
We leave the detail expressions of these WZNW fields in Appendix A. With the help of 
Eq. (|A1|) and shifting $ s . — ► $ s + 2ao X i ^ ne correlation functions in the presence of a weak 
magnetic field can be shown as follows: 



(2?(r)27(0)> = (77(t)T;(0)) = (2?(r)T/(0)> 

1 1 



+ 



47T 2 (VfT + ix) 2 (VfT — ix) 2 

{cos + M)x] + cos [£(1 - M)x]} 



+(-l)«o 



M 2 



( v 2 r 2 +x 2)l/4 a ( t; 2 r2+a; 2)l/2 
A 2 



V 2 T 2 + X 2 ' 



1 



+ 



1 



{cos + M )x] + cos [£(1 - M)x]} 



2ap 



]t 2 + x 2 y/ 4a (v 2 r 2 + X 2 ) 1 / 2 



+A 3 



cos [f Q (l- M)x\ 



{v 2 r 2 + x 2 y/ a 



(St(T)S7(0)) 



(V 2 T 2 + X 2 )^/ 2 



i^M x VsT -i x \ _ i ^ Lx / VsT + j x 

r — ix 



VfT + IX t 



+ 



{V 2 T 2 + X 2 yl\v 2 T 2 + X 2 ) 1 / 2 



+ (-l)«0 



{ V 2 T 2 + X 2 )^/ 2 



iZM x / Vs T - IX \ ^ILM X / V S T + IX 

g 2a 



_|_ g 2a ' 



y fT — ix t 



VfT + IX t 



(10) 



where r is the imaginary time and x = a \i 71 = 1 + | + A- > 2. Aj and -Bj (i = 1, 2, 3) 
are nonuniversal constants. 

The various characteristic momenta in Eq. ( |T0"D can be understood as the combination of 



S 



the shifted Fermi momenta under the magnetic field: fepf = ^(l+M) and kp\ = i^(l — M). 
The structures of orbital (spin) correlation functions are not affected by the magnetic field 
(uniaxial pressure) except the 2kp component. The original degenerate 2/cp-excitations now 
split into two incommensurate soft modes at 2k p] = <^(1 + M) and 2k f\, = <^(1 — M). 
Moreover, the corresponding exponent r )2k F = 1 + > 1-5 even at M = + and J = K/4. 
{l2k F — 1-5 when M = and J = K/4.) This is because the marginally irrelevant coupling 
G3 in Eq. @ now turns into a marginal one and thus changes the compactification radius 
of $ s as shown in Eq. (§). The 4fcp^ mode does not appear in (T i a (r)Tj 3 (0)). (Note 
that 4fcpj = 4kp j [ (mod. 2tt) and thus the 4kp^ mode is not independent.) This means 
that it does not carry nontrivial orbital quantum numbers. The corresponding exponent 
7 4 fc F = 2/ a = 4(7 2 fc F — 1), which is consistent with the prediction of conformal field theory 
[ 11 1 . The only enhanced fluctuation is the 2kp component of [Sf{r) Sj(fy) (the B 2 term) 
of which the exponent becomes 1 + a/2 < 1.5. Another interesting observation is that the 
position of the cusp at q = tt/2 in the static transverse spin structure factor does not shift 
under the magnetic field. 

B. J > K/4 

When G > 0, there are two massive regions. As suggested in Ref. @, the Hamiltonian 
(H) with J > K/4 falls into the region with < G < —2Gs in the absence of the magnetic 
field. We expect that under a weak field it is still the case which corresponds to the region 
II. In this region, because g + is marginally irrelevant and the 0(4) symmetry is restored as 
g + = (G = —2G 3 ), we conjecture that up to logarithmic corrections, H so is equivalent 
to the 0(4) GN model at strong coupling regime. Next, we examine the effect of H int on 
H s . Integrating out £ a (a = 3, ■ • ■ , 6) gives the interaction term — ^ cos (\/4jt& s ) where 
m = K92CEI=3^l) + 2 #i (£!££))• ( Tli e term (<9 Z $ S ) 2 is possibly generated. But its effect 
can be absorbed into the g\ term in H s .) By replacing H s with the following one: 

H s = fdx \ V ^i(d x e s ) 2 + (d x $ s ) 2 }-— cos(v / 4T<I» s ) + ^(^<I> s ) 2 - \ (ii) 
J [2 rcao ix \Jix J 
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where v s has been substituted into to incorporate its renormalization effect, we assume 
that the low energy and long distance behavior of the Hamiltonian ([I]) with J > K/4 and 
h > h c is approximately described by H so and H s . 

By introducing a Dirac fermion i/)r(l), H s can be fermionized as: 



H s = dx [-iv s {ip^d x tp R - i)td x ip L ) - im(^ij; L - ipt^n) ~ K^r^R + ^1^- 



Lj 



+ 9 -±^ r + ^l)\ (12) 

7T 

We see that h is equivalent to the chemical potential of fermions and h = h c = \m\ is a 
quantum critical point. We are concerned with the case where h > h c . In this case, the low 



energy excitations of the $ s sector are described by the following effective Hamiltonian pi 

1 



H — dx — 
2 



~{d x 6) 2 + g{d x <j>) 2 
9 



(13) 



where is obtained via shifting $ s by an amount 1 x. g is a parameter determined by gi 
and g < 1 (g > 1) when g± < (gi > 0). The scattering (gi) term has only negligible effects 
in the limit M — > where we will get free fermions with g — > 1. 

Based on our assumption that H so is equivalent to an 0(4) GN model in the strong 
coupling regime, we replace H so with the following one: 

Hso = -- v f - + g(t ^) 2 ' ( 14 ) 

a=3 a=3 

where ^ > 0. We have performed the duality transformation: £ R — > £ R , ^£ — > — ^| , and <r 6 <-> 
/ig. The same transformation must be applied to Eq. (|Al]) before calculating correlation 
functions. The spectrum of 0(4) GN model is distinguished from that of the O(N) GN model 
with N > Q |[22j| . The latter consists of the elementary fermions, the bound states of them, 
and kinks and can be qualitatively captured by the large N approximation. The existence 
of kinks is intimately related to the spontaneous breaking of discrete chiral symmetry (Z 2 
symmetry), which corresponds to the breaking of lattice translation symmetry in our case. 
On the contrary, in the 0(4) case, the elementary fermions are unstable against decay 
into kinks which become the only stable excitations. As a consequence, the magnetic field 
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dramatically changes the spectrum of the S z = sector. As shown in Ref. ||, the low energy 
excitations in the absence of the magnetic field are massive fermions which belong to the 
representations (Stat,T tot ) = (1,0), (0,1). Under a magnetic field, the symmetry becomes 
U(l)xSU(2) and the good quantum numbers are S z and T tot . As we expect, the excitations 
of $ s carry (S z ,T tot ) = (±1,0). Nevertheless, those massive fermions carrying S z = are 
no longer stable when h > h c . The elementary excitations turns out to be kinks which carry 
(£", I 1 **) = (0,1/2). 

A convenient way to deal with Eq. ([14]) is to transform it into two decoupled sine-Gordon 



models [23 



H so=J2{lT j(d x e a ) 2 + l(d x ^ a ) 2 -44cosVM a L (15) 



a=± 



2 



1 
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7r 2 ag 



where I = Jl + anc [ ^ — Since I = 1 + , the coupling constant of the sine-Gordon 
model (3 2 = 8tt~ . The kinks of the 0(4) GN model are nothing but the solitons of the two 
sine-Gordon models. Within our approximation Eqs. (|13|) and (0), we are able to discuss 
the dynamical structure factors of spin and orbital sectors, which are defined as: 

/oo 
dxdr (S«{T)Sf(0))e iuJT - iqx , 
-oo 

/oo 
dxdr {T^{r)Tf{Q))e iu>T - iqx , (16) 
-oo 

where x = ao\i — j\. For T a/3 , it suffices to compute T zz because the orbital SU(2) symmetry 
remains intact. We shall see that the 2/c^ components depend on x + y as well as x — y. 
This is a manifestation of the spontaneous Z2 symmetry breaking of the ground state. In 
that case, we just list the correlators in the coordinate space. The details of computations 
are left in Appendix B and we show the results in the following. First, we consider the spin 
correlators. They are 

l— 1 



S zz {u, q) = -z- u[5{uj - v s q) + 6{u + v s q)} + C x 



gv s ' ' s \uj 2 - v 2 s {q - Ak Fi ) 2 j 

x[Q(uj - v s (q - 4k Fl ))Q(u + v s {q - Akp^)) - (w -> -u)], 



Wr W »))~™(^)«(i,)(g) , (17) 
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where x = iao, y = jcto, and z = v s r + i(x — y). C\ is a nonuniversal constant. For S zz , 
the C\ term is correct only near its low energy threshold and q ~ 4fcp^. In addition, it 
diverges at the low energy threshold as [uj ±v s (q — Akpi)Y 1+ '9 for g > 1 while it approaches 
zero with the same functional form for g < 1. As M ^ 0, the exponent becomes zero. For 
(S^~(t)SJ' (0)) , as the case with J < K/4, the characteristic momentum q = n/2 does not 
shift under the magnetic field. Moreover, the exponent approaches ~ as M — > 0. 
Next, we consider the orbital correlators. The results follows: 

T (w ) ?«0) = 8a ' sgn(u), 

s J vs — 4m ^ 

(^(r)T/(0)}| 2fcF ~ 2 cos [£-(x + y)] cos [^(s - y)) 
x[Wi(T,x-y)-W 2 (r,x-y)] 
+{cos [2k F ^(x — y)\ + cos [2kp^(x — y)]} f •py 



j_ 



x [Wi(t, x -y) + W 2 (r,x - y)}, 

T**(u, g « * ) = C 2 , ^ x sgn(u), (18) 

V a W« 2 -4m 2 cos {f) + cosh [fx(w)] 

where x = ia , y = jao, 2; = 5 s r + i(x — y) and m is the mass of kinks, s 2 = uj 2 — v 2 q 2 , 
u 2 = u 2 - vj(q - n/a ) 2 and C = Trf3 2 /(8ir - f3 2 ) with (3 2 = 87T". x(s) = cosh_1 (s/2m). 
The functions /(x) and F (x) are defined in Eqs. (|B7|) and ( |B11|) , respectively. C 2 is a 



mere constant. The result for T zz (u,q ~ 0) is exact when 4m 2 < s 2 < 16m 2 while that 
for T zz (u,q ph ^) is valid when 4m 2 < u 2 < 16m 2 . For higher energy, there are small 
contributions from four, six, eight, etc. particles states. The functions Wi }2 , which are 
defined in the following: 



Wi(t, x) = F(sin h/7r/2($+ + $-)] sin Wn/2(Q+ - 6_)]) 



W 2 (t, x) = F(cos [ v/7r/2($ + + $_)] cos [yJir/2{Q+ - 6_)]), (19) 

with F(0) = (0(r, x)O(0, 0)), involve the computation of form factors containing <&± and 
their dual fields Q±. At present, no one knows how to compute them. However, we can still 
discuss some features of these functions. The operators we are considering are like A + ■ 

12 



where A± = cos [or sin](W7r/2$±) x cos [or sin](W7r/26±). According to Ref. P7|, the 



operators cos [or sin](y 7t/2Q±) are fermionic while cos [or sin](y ir/2Q±) are bosonic for 
(3 2 = 8tt~ . Thus, A± must be fermonic. This implies that the first nontrivial form factors of 
A± start from the one-soliton states. As a result, the leading contributions to Wi come from 
two-particle states. We expect that in the momentum space, the 2kp component exhibits 
incoherent background with two-particle thresholds. 

In summary, when h > h c , the spin correlators exhibit algebraic decay while orbital 
correlators remain massive behavior. Low energy modes appear close to q = and 4k for 
gzz QT g = jl_ f or 5*+- 2k p component of S zz and the uniform and 4kp components 

of S + ~ involve massive excitations [see Eq. ("Kl)] and thus exhibit similar behavior to 



that of (T?(t)T?(0)) \ 2 k F - The leading contributions to the uniform and 4kp components of 
(T"(t)T?(0)) come from the two-particle states, which is similar to that in the absence of the 
magnetic field. The main distinction between the cases with h > h c and h < h c is the 2k p 
component. For the former case, it also exhibits two-particle thresholds in contrast to the 
latter case where it shows a coherent "magnon" peak. This can be detected experimentally 
by examining the dynamical spin structure factors under a uniaxial pressure. 

IV. WEAK COUPLING REGIME: K < J 

When K J, a proper starting point is to consider the K term in Eq. (fj) as a 



perturbation, which was discussed by Nersesyan and Tsvelik [16] in the context of two- leg 
spin ladders. The unperturbed Hamiltonian is composed of two antiferromagnetic (AF) 
Heisenberg chains which can be treated by using the standard relation between the spin 



(pseudospin) operators and WZNW fields P0| : S n /a = j s + (-l)"om s and T n /a = j t + 



(— 1) a o fh t where x = na^. By defining <p± = (<f) s ± (f> t ) / V2 where (p s and <p t are, respectively, 
the bosonic fields describing the spin and orbital degrees of freedom, we obtain 



#=E fdx {li(dJ a ) 2 + (dM 2 } + —cos(^ a )-^=d x 
^ ± J [2 ira V47T 



, xfa ( j 



(20) 
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where v ~ Jciq and m = X 2 K/(2ir 3 ). The Hamiltonian ( p0|) can be further fermionized as: 



-^aR^aR + ^tL^PaL)}- 



(21) 



From Eq. (|2lD , we see that in the region with K <C J, the low energy excitations in 
the absence of the magnetic field are described by two free massive Dirac fermions or four 
massive Majorana fermions, which are kinks connecting two degenerate ground states ||16|| . 
In our case, these Majorana fermions belong to the representation (S to t,Ttot) = (1/2, 1/2) 
and are coupled to the magnetic field simultaneously. This leads to the result that the 
QCP moves to h c = 2m instead of m. Beyond that value, all excitations become gapless 
and the correlation functions are algebraic decay. The situation is quite different from the 
two-leg spin ladder model by considering Tj in Eq. ([[]) as the spin operator on the other 
chain. In that case, the magnetic field is only coupled to ip + and the QCP is located at 
h c = m. Beyond this value, ip- is still massive and we expect that the behavior of correlation 
functions is similar to that of the usual spin ladders. 

What we are concerned is the case where h > h c . The effective Hamiltonian describing 
the low energy excitations around the Fermi point is as the following: 



a=± J * 

Note that Eq. (^) is in fact a theory describing two free massless fermions. We are now in 
a position to calculate the spin and orbital correlation functions. In terms of <j) a and 6 a , the 




(22) 



spin and orbital operators are expressed as follows PTf : 



- + -^,(0 + + 0„) + A 1 (-irsin[ 

L A/7T 



irM 



x + v / 7r(0+ + 0-)], 





+X 2 (-l) n e i ^ i§++§ -\ 



rpz 

n 



2^ 



d x (<p + - 0_) + A 3 (-l) n sin [v^F(0 + - 
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r+ = i- 

Ztt 



+A 3 (-l) n e^+-^-), (23) 

where x = na^ and Aj (i = 1, 2, 3) are constants. Note that S* 2 contains the Q (= ^-(1 — M)) 
term in contrast to the two-leg spin ladder where the term next to the uniform component 
is the 2Q one [21]. The reason is that in the spin ladder, the antisymmetrical mode 0_ is 
still massive even under the magnetic field and leads to the exponential decay of the Q term. 
However, in our case, 0_ becomes massless when h > h c . With the help of Eq. (|23|) , we get 
the correlators as the following: 

WWS?(0)> = ^ + ^ (i + + B,«.(Ox)^, 

o _ 

(£+(r)Sr(0)) = — 5_ U~ x - + e-' l - x -) + 



z ZJ \z\ 



(2?(r)27(0)> = (^(r)T/(0)> = (^(r)^(0)) 



?7T 2 VZ 2 Z 2 / IZ 



where x = ao\i— j\, z = vr+ix and Di (i = 1, 2, 3) are constants. Note that all correlators in 
Eq. (pD behave like the spin correlation functions of an AF Heisenberg chain when h — > h+ 



(M — > + ). Low energy modes appear close to g = and Q for S 122 or g = ^ and ^ for 



S + while for spin ladders they are close to q = and ^ for S zz or q = Q and for 



S + The orbital correlators are identical to the spin correlation functions of an AF 



Heisenberg spin chain up to some numerical prefactors for any finite h (> h c ). 



V. CONCLUSION 

In this paper, we discuss the effects of a magnetic field on the one dimensional spin- 
orbital model with J\ = J 2 = J. In the gapless phase, i.e. J < K/4, the spin and orbital 
correlation functions still exhibit power-law decay. The magnetic field manifests itself on 
the emergence of incommensurate soft modes with characteristic momenta depending on 
the magnetization. Furthermore, the exponents are also dependent of the magnetization 
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though they have universal values at zero field. A distinctive feature in the case with orbital 
degeneracy is that the structures of orbital correlators are hardly influenced by the magnetic 
field. This is associated with the unbroken orbital SU(2) symmetry. 

In the massive phase, i.e. J > K/4, we expect that there is a quantum phase transition 
induced by increasing the magnetic field. However, because of the different excitation spec- 
trum at different range of J / K, the correlation functions exhibit distinct behavior between 
the region near the SU(4) point and the one with K <C J. In the former case, the magnetic 
field exerts a great influence on the S z = sector such that the spectrum is completely 
different from that in the absence of the magnetic field. We propose that this change can be 
experimentally examined by studying the dynamical spin structure factors under a uniaxial 
pressure. In the latter case, the spin and orbital correlation functions become power-law 
decay when the magnetic field is greater than the critical value. Furthermore, the exponents 
take universal values which are identical to those of an AF Heisenberg chain. This is due to 
the fact that the underlying system is described by a free fermion theory. 

Finally, we would like to mention the applicability of our results to the asymmetric case 
( Ji 7^ J-z) around the SU(4) point. The massless and massive phases on the symmetric line 
now extend to large anisotropic regions JI3| . The phase boundary passes through the SU(4) 
point. The anisotropy affects not only the coupling constants in the Hamiltonian (|7|) but 
also the velocities of the spin and orbital sectors such that they become different. For the 
massless phase with G 2 = c(J 2 — K/A) < (c > 0), the correlators under the magnetic 
field are similar to Eq. ( fLOl ) except that the velocity of £ 6 is different from the one of the 
orbital sector. The symmetry of the fixed point Hamiltonian is in general U(l) xSO(3) xZ 2 
instead of U(l)xSO(4). In the massive phase with G 2 > 0, we still have an approximate 
U(l)xSO(4) symmetry at low energy for h > h c and the correlators are the same as Eqs. 
(|T7|) and fllSp . On the other hand, in the rest part of the phase diagram around the SU(4) 
point, it is possible that another phase transition of the KT type occurs on the S z = 
sector when h further increases However, it is just a transition from the massive (our 
region II) to massless (our region I) phases. Therefore, the correlators are still of the forms 
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we obtained in the corresponding regions. 



ACKNOWLEDGMENTS 

Y.L. Lee is grateful to P. Lecheminant for enlightening discussions. The work of Y.L. 
Lee and Y.W. Lee is supported by National Science Council of R.O.C. under the Grant No. 
NSC89-2811-M-007-0015. 

APPENDIX A: THE SPIN AND ORBITAL OPERATORS 

In this Appendix, we list the expressions of the WZNW fields J a (t), N s (t), and n s (t) in the 
following: 

% - 1 

Jtu — ~ itu x £tu, Jl — —^=d x ^ s , 
2 WTT 



J s+ = --^L fcexp [-i^(<S> s - Q 8 )] + ^exp [i^s + Os)}} , 

N* = A[sm (V7r$ s )A* 3 A*4^5A*6 + i cos (v^^ow^oe], 
N s ± = v477i773(/i3/i 4y U50- 6 ± 0-30-40-5/16) exp (±iy/7rQ s ), 

N* = A[cos (^^ s )a 3 a 4 n 5 fx e + i sin WaWb], 
N t ± = -iAr] 2 r] S (a 3 iJ,4 ± 1^304) [sin (v^*)/^ =F cos (V^S,)^/^], 

n t = -ifiSfl x £ tL , = sin (v 7 !^), 

7rao 

n s+ = ^ {e| exp [iV^(^ s + © s )]+^i exp [-iV^(* s -©,)]}, (Al) 

where u = R,L and 0± = O 1 ± iO v (O = J s (t), N g (t), n s (t))- A and B are nonuniversal 
constants. £ tw = (z£^, — is a vector under orbital SO(3) rotations. <r a and /i a (a = 

3, • • • , 6) are order and disorder parameters of the corresponding Ising models. 

APPENDIX B: CORRELATION FUNCTIONS IN THE REGION WITH J > K/4 

Following the appendix of Ref. |f2TJ, we found that in this case it is not necessary to 



include higher harmonics. When h > h c , the spin operators become 
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M an M , 

St = — + ^=d x <j> + A(-l)«o sin (—x + v^0), 

I y/TT a 



~ e'^ 1 ^ cos (v / 2^$ T ) + H.c, (Bl) 
where A is a constant. As for the orbital operators, the corresponding WZNW fields become 

y In 

+i sin (— x + v^) cos - ($+ + $_)] cos - (0+ - 6_)], 
n t z ~ sin [v / 2^($+ + $_)]. (B2) 

The computations involving the <ft sector are straightforward. To compute the longitudinal 
dynamical spin structure factor, we need the following integral: 

e iu)t-iqx 2-ji- 2 e i7T ( a ~P)/ 2 

-oo ^ (x + vt-iO+nx-vt + iO+y = @( " - VqMu + n) vT(a)m 

(B3) 



2v V- a { 2v 



\uj — v q J \io + vqj 
The correlators involving $± only can be calculated in terms of the exact results of form 
factors on the sine-Gordon model [p4]-|27l]. The calculations can be partially simplified by 
noting that the Hamiltonian fll5|) is invariant under the exchange $ + «-> $_ and $ + and $_ 
are decoupled. We now compute T zz {oj, q ~ 0). (A similar calculation has been done in the 



context of a spin ladder model |29fl .) With the help of Eq. (p2|), it is 

T zz (lu, q w 0) = — ^ Im lim / dzdr (9 a; $ + (r, x)a a; $ + (0, 0))e^ r -^ 

/OO /'OO . 

dx ^e^ +j0+ )*-^([jo(t,x),jo(0,0)]), (B4) 
-oo JO 

where jo = 77^^+ is the temporal component of the current operator in the sine-Gordon 
model and we set Vf — 1. In one dimension, energy and momentum can be parametrized 
by rapadity 9 as e = m cosh 9 and p = m sinh 9 where m is the mass of the sine-Gordon 
solitons. The resolution of the identity is given by 

oo . jn jq 

i = EE / ,'\' n , \ n, o 1 ) an ... ai a —(9 1 , ...,e n \, (B5) 



n=0 ct. 



(2vr) n n! 



where n is the number of particles, = ~, — \ for solitons and antisolitons, respectively. 
Inserting Eq. ( |B5|) into Eq. ( p34|) and performing the integrations over space and time yields 



r>, q « o) = -4™g im EE/ ^z^l^ ,.^^, • • • A)l 2 

n=0 Oj J ' t - 



X 



i(?-m^; sinh 



S(q + mJ2j sinh 6* j 



uj — m J2j cosh 9j + i0 + u + mJ2j c °sh Qj + ^O 4 " 



where 



f 3 ° (9i, ■ • ■ A) = (o|jo(o, o)\e n , • ■ • , ^} ari ... Q1 



(B6) 



is the sine-Gordon current form factor. (For an operator d(t,x), the form factor 
F °,...aJPu'''M = (O|O(O,O)|0 n , ■■■,9i) a n —ai-) Note that an n-particle state only con- 
tributes to Eq. ( |B~6| ) when s 2 = uj 2 — q 2 > n 2 m 2 and n must be an even integer. Thus, at 
low energy s 2 < 16m 2 , only two-particle states contribute. The corresponding form factor 
I ?3 is 

9 1 + 9 2 * 



2 • 2 



-2m sinh 
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where 6 I 12 = #1 — 9 2 and 

/(f) = i~ s^h (0/2) x exp { f°° *L ^ [{ ° " .^ /2] [tanh (ttx/2) - 1]). (B7) 
Z7r [Jo x smh (7TXJ J 



After performing the integrations over 9i, we obtain 

r-( M ,, w o) = ^J^«»'V( M ), 

s^v s z — 4?7r 

where x(s) = cosh -1 (s/2m) and 4m 2 < s 2 < 16m 2 . 



(B8) 



Next, we compute T zz (oo,q rs ^). According to Ref. p8| , (cos(|0)) = Qp/2 where 



mVir(47r/(87r - j3 2 )) 
2r(/3 2 /(167r - 2/3 2 )) 

xexp U T 



4tt 



sinh 2 (2a(3t) 



-16nt 



2 sinh (f3 2 t) sinh (8vrt) cosh ((8tt - (3 2 )t) An 



(B9) 
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and P < \^8n is the coupling constant of the sine-Gordon model. Therefore, the leading 
term of (nf (r, x)n z t (0, 0)) oc Qj /2 F{sm (f $+)) where /3 2 = 87T - . Then, 

T zz (uj,q^ — ) oclmz /°° dx /°° dte <(a,+i0+) *" i(5 "^ )a: ([sin $+)(*, x), sin $+)(0, 0)]) 



S(q — mJ2j sinhfy 



5(q + mJ2j sinh 9j 



to 



(BIO) 



The leading contributions to Eq. ( BIO ) come from the two-particle states and the corre- 
sponding form factor [27J is 



sin(f$+), 



12 



7 12j 



'l2J OC 



W12J 



COsh [|?(012 - «7r) 



where ( = vt/? 2 /(8tt - /? 2 ) and 



F (6) = -i sinh (9/2) x exp 



00 da; sinh 



•7T-C 



;r 



sin 2 [f(6> -«tt)] 



o x sinh (7rx) sinh (^) cosh (^) 



(Bll) 



Here the indices of the form factor a, = 1,2 correspond to the neutral fermions which are 
related to the sine-Gordon solitons through Z± = Z\ ± %Z 2 where Z± and Z 12 are, respec- 



tively, annihilation operators of solitons (anti-solitons) and neutral fermions |27|. Inserting 
Eq. (pil|) into Eq. (|B 1Q ) and performing the integrations over 9{, we get 

\F (2 X (u))\ 2 



T (u, q « — ) oc — ^=^= x ; 

ao MV« 2 - 4m 2 cos (^-) + cosh [*x(«)] 



sgn(u) 



(B12) 



where 4m 2 < u 2 = uj 1 — (o — 7r/ao) 2 < 16m 2 . 
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